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Abstract 

We apply the Born-Oppenheimer approximation to a harmonic diatomic molecule with one 
electron. We compare the exact and approximate results not only for the internal degrees of 
freedom but also for the motion of the center of mass. We address the problem of identical nuclei 
and discuss other applications of the model and its limitations. 
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I. INTRODUCTION 



The first step in the treatment of a classical or quantum-mechanical isolated system of 
particles should be the separation of the motion of the center of mass from the internal 
degrees of freedom. Most textbooks on quantum mechanics and quantum chemistry do that 
for the hydrogen atom but then consider the nucleus at rest when they discuss many-electron 
atoms 1,2 . Such omission is more marked in the case of molecules because they invariably 
resort to the Born-Oppenheimer (or clamped-nuclei) approximation (BOA) and do not even 
mention the problem at alR However, an adequate treatment of the motion of the center of 
mass is most important for the estimation of adiabatic corrections to the potential-energy 
hypersurface 3 

In two earlier discussions of the BOA in this journal the authors chose two-particle mod- 
els: the hydrogen atorn^ and a delta potential in a box-. However, the simplest molecule, the 
hydrogen molecular ion Kj~, is a three-particle system. Therefore, those earlier discussions 
are not sufficiently realistic to show many relevant features of the treatment of molecular 
systems and of the nature of the BOA. 

The purpose of this paper is the discussion of the separation of the motion of the center of 
mass in the BOA. As already mentioned above, this aspect of the problem is not discussed 
in most textbooks on quantum chemistry- and the BOA is almost entirely omitted from 
those on quantum mechanics^. For simplicity we resort to a simple model of three particles 
(two nuclei and one electron) that interact through Hooke's forces. Although such harmonic 
interaction is unphysical, it has the advantage that the Schrodinger equation is solvable 
analytically and we can thus compare the exact solution with the exact BO one. 

In Sec. [Til we present the model and obtain a dimensionless Hamiltonian operator in a set 
of variables that allows the exact separation of the motion of the center of mass. In Sec. IIHI 
we solve the Schrodinger equation exactly. In Sec. [IV] we obtain the exact BO solution and 
compare it with an expansion of the exact nonadiabatic result derived in Sec. IIHI In Sec. [V] 
we address the interesting case of identical nuclei and discuss the permutational symmetry 
of the molecular states and the correlation functions. Finally, in Sec. I VI I we discuss other 
applications of the harmonic model and its limitations. 
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II. MODEL 



Our model consists of three particles of masses Ml, M 2 and M 3 that move in one dimen- 
sion and interact through forces that follow Hooke's law 

2M x dXl 2M 2 dX$ 2M 3 dXj ^ ' 



+5 



K a (X 1 - X 2 f + K l3 (X 1 - X 3 f + K 23 (X 2 - X 3 f (2) 



where are the force constants. 

In order to make the Schrodinger equation dimensionless we choose units of length L, mass 
M and force constant K and define the corresponding dimensionless quantities Xi = Xi/L, 
mi = Mi/M and kij = K^jK. If we choose L = h 2 /(MK) ^ and take into account 
that h 2 / (ML 2 ) = hcu, where cu = \Jk/M, we easily derive the dimensionless Hamiltonian 
operator 

H Id 2 Id 2 Id 2 



huj 2mi dx 2 2m 2 dx\ 2m 3 dx^ 



+ 2 



k u (xi - x 2 ) + fci 3 (xi - x 3 ) + k 23 (x 2 - z 3 ) (3) 



From now on we write H instead of H /(Tiu) and simply remember that the energy is mea- 
sured in units of Hcu. Notice that we may choose M — Mi and/or K = in order to 
make some particular parameters equal to unity and simplify the Schrodinger equation even 
further. However, we have decided to keep all the model parameters for clarity. 

In order to separate the motion of the center of mass from the internal degrees of freedom 
we define new coordinates as follows: 

qi = — (niiii + m 2 x 2 + m 3 x 3 ) , m T = m 1 + m 2 + m 3 
rriT 

q 2 = x 2 - xi 

q 3 = x 3 - xi (4) 

Notice that qi is the coordinate of the center of mass and q 2 and q 3 are simply the positions 
of the particles 2 and 3 with respect to particle 1. The latter variables are translationally 
invariant as it is expected for the internal (spectroscopic) degrees of freedom. More precisely, 
the displacement of the whole system Xi — > Xi + a produces q\ — > qi + a, q 2 — > q 2 and q 3 — > q 3 . 
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The Hamiltonian operator ([3]) (in units of hu) becomes 

Id 2 Id 2 1 



H 



d 2 



2rriT dq\ 
1 



2m 2 dq 2 



2m 3 dql 



d 2 



2mi 
1 



d 2 d 2 

I |- 2 

dq\ dql dq 2 dq 3 



+ - 



(5) 



III. EXACT SOLUTION 



The Hamiltonian operator (J5]) is the sum of an operator that depends only on qi 

H, = _ J_^. 

2m T dq 2 

and another one that depends only on q 2 and g 3 



(6) 



Ho 



1 d 2 



1 «9 2 



1 



d 2 



d 2 d 2 
+ 2-, 



2m 2 dq 2 2m 3 <9g§ 2m\ \dq[ dq% dq 2 dq, 



+ fc 12?2 + k 13ql + h3 (?2 - 93; 



(7) 



Therefore, each eigenfunction of H can be written as the product of an eigenfunction of Hi 
times an eigenfunction of H 2 . The operator Hi describes the motion of a free pseudoparticle 
with mass equal to the total mass of the system. If Q(q 2 , g 3 ) is an eigenfunction of H 2 



we conclude that the eigenfunctions ^(91,92,93) of H 



Hip = Etp 



(8) 



are of the form 



and 



^(91,92,93) = e ift<?1 $(g 2 ,g 3 ) 



E 



2m<r 



+ e 



(9) 



(10) 



(11) 



where —00 < k < 00. 

Notice that the Hamiltonian operator H has a continuous spectrum, and that the Hamil- 
tonian for the internal degrees of freedom H 2 has a discrete (or point) one (which we will 



obtain below). This point is most important when one has to calculate expectation values 
of observables. For example, everybody knows that if the potential-energy function is ho- 
mogeneous of degree two (that is to say, it satisfies V(tr) = t 2 V(r)) then the virial theorem 
gives us<T>=<y> for an eigenfunction of the Hamiltonian operator. However, this 
relationship does not apply to the kinetic and potential energies in H because the integrals 
diverge, but it applies to the kinetic and potential energies in H2 since its eigenfunctions are 
square integrable. 

The Hamiltonian operator ([Tj) is a particular case of 

* = -S^4 + 5??^ (12) 

In order to obtain its eigenfunctions and eigenvalues we carry out a change of variables of 
the form 

?< = 2cii%- ( 13 ) 
3 

Thus the Hamiltonian operator ffl2l) becomes 



1 



T 



d 2 1 



h = -^EE c-a(c-) ..^ + ,EE(c r BcL^. 



where A, B, and C are matrices with elements Aij, B^, and Cy, respectively. We choose 
the matrix C in such a way that 

c^Ajcr 1 ) 7 = I 

C T BC = D (15) 

where D is a diagonal matrix 

Dy = UJ 2 5ij (16) 

This approach is well known in the treatment of small oscillations in classical mechanics 6 . 

The resulting Hamiltonian operator is a sum of uncoupled dimensionless harmonic oscil- 
lators 

d = "5?S + 5?"M (17) 

therefore its eigenfunctions are products 

${n} = II^(2/i) ( 18 ) 
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and its eigenvalues are given by 

€ {n} = H ' e «* = ( n i + 2) Ui ^ 
where rii = 0,1,... are harmonic-oscillator quantum numbers, and each a;, is a dimensionless 
frequency. More precisely, 

(id 2 1 \ 

[~2dy 2 + 2^ y V = en ^nAVi) (2°) 

It follows from equation ffl5|) that 

C X ABC = D (21) 

which shows that the problem reduces to the diagonalization of the nonsymmetric matrix 
AB. 

In the particular case of the Hamiltonian operator ([7j) we have 



A 



B 



/ j_ 1 j_ j_ 

mi rri2 mi 

— — + — 

\ mi mi m3 

' &12 + k 23 -^23 ^ 

-^23 fcl2 + ^23 



fci2(mi+m2)+A,'23mi fci 3 m2— &23ra 



AB = mim2 mim2 (22) 

! k 12 m s -k 2 : i m 1 fci3(mi+m 3 )+fc 2 3mi 1 

The characteristic polynomial for the matrix AB is 

2 k 1 2 , m 3 (m 1 + m 2 ) + k 1 3'm 2 (m 1 + m 3 ) + /c 23 mi(m 2 + m 3 ) 

w — w 

mim 2 m 3 

| (mi + m 2 + m 3 ) [fci 2 (fci3 + ^23) + ^13^23] ^3) 
m\m 2 m 3 

where w = uo 2 . The two real positive roots give us the frequencies that we need to obtain 
the energy eigenvalues according to equation (1191 : 

€n u n 2 = Wl (ni + -^j+UJ 2 (n 2 + (24) 

We do not show those roots explicitly here because they are rather cumbersome. 
The exact eigenfunctions for the internal degrees of freedom 

^ni ,n 2 (92, Q3) = W (2/2)^2 (1/3) (25) 

clearly show the coupling of the motion of the particles through the variables jji that are 
linear combinations of the qj. In other words, the problem is completely separable in the 
variables but not in the qj or 



IV. THE BORN OPPENHEIMER APPROXIMATION 



The BOA is discussed in many textbooks 2 and also in earlier pedagogical articles in this 
same journal 4,5 . For this reason we do not develop it here explicitly and just show its results 
for the present model. From now on we assume that our three-particle system models a 
diatomic molecule with just one electron. We choose the particles 1 and 2 to be the nuclei 
and particle 3 to be the electron; more precisely, we assume that m\ > m 2 ^> m 3 . We clearly 
appreciate that the harmonic potential chosen here is not realistic because it describes an 
attractive interaction between the nuclei. However, since the BOA is based on the different 
particle masses and not on the nature of the interaction we can apply it successfully and 
compare its approximate solutions with the exact ones. 

In the clamped-nuclei approximation we omit the kinetic energy of the nuclei, which we 
assume to be at rest at x\ and x 2 , and solve the Schrodinger equation for the remaining 
"electronic" Hamiltonian 2 

1 d 2 



2m 3 dx 3 



+ 2 



&13 (Xl ~ X 3 ) + k 23 (x 2 - X 3 y 



(26) 



Since the "nucleus-nucleus" interaction k\ 2 {x\ — x 2 ) /2 is just a constant we add it later 
to the electronic eigenvalues. 

If we rewrite the potential-energy function as 

l 2 



ki3 (xi - x 3 ) + k 23 (x 2 - x 3 
1 



(fcl3 + k 



2?,) 



%3 



k i3 xi + k 23 x 2 



k 13 + k 



28 



ki 3 x ± + k 23 x 2 



(kyjXi + k 23 x 2 f 
2 (k 13 + k 23 ) ^ 2 

then we realize that the electronic Hamiltonian (1261) is just a displaced harmonic oscillator 
and that the electronic energies are given by 

- ^ + ^ + \ (k 13 x\ + k 23 xf 
2 (k 13 + k 23 ) 2 V i 23 2, 



e e,m X 2 ) 



k 13 + k 



m 3 



rii + 



(27) 



The nuclear motion is governed by the potential-energy function 

1 2 

U(xi,x 2 ) = U(xi - x 2 ) = e e , ni (xi,x 2 ) + -kn (xi - x 2 ) 



k 13 + h 



2\\ 



m 3 



ni + 



. k\ 3 k 23 + k\ 2 ki 3 + k\ 2 k 23 2 
H — — ; — ^ (Xi - x 2 ) 



2 (fci 3 + h 



23 I 



(28) 
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The final step is the solution of the Schrodinger equation for the nuclear Hamiltonian oper- 
ator 



Hn = -t^t^ - jr-jr* + U{x x - x 2 ) (29) 



1 d 2 

2m\ dx\ 2m 2 dx\ 

In doing so, we separate the motion of the center of mass by means of the change of variables 



R = (rriiXi + 1712X2) , rriN = mi + W2 

m N 

q 2 = x 2 - X! (30) 

and rewrite the Hamiltonian operator (1291) as 

Tj Id 2 Id 2 mim 2 

Hn = ~^r N dR 2 ~ 2ir N M + u(q2) ' ^ N = ^r (31) 

The eigenf unctions of this operator are of the form 

<S> N {R A2 ) = e iKR y{q 2 ) (32) 
and its eigenvalues are approximations to those of H: 



F BO ^ /fcl3+_^23 / 1\ ^13^23 + ^12^13 + kukq 

K ™ ~ 2m N + \ m 3 r 1 + 2 ) + ]] (k 13 + k 23 ) fi N V" ' 2 

If & e (x3] Xi, X2) is an eigenfunction of the electronic Hamiltonian H e then the full BOA 
eigenfunction is 

V^Gzi,?^) = e iKR ^ e (x 3 ;x 1 ,x 2 )< f (q 2 ) (34) 

Although the BOA treats the electronic and nuclear motions separately, we appreciate that 
the BOA solution (1341) already couples the nuclear and electronic degrees of freedom through 
the electronic function ^ e {.x^ x\, x 2 ) much in the way the exact solution (1251) does. Besides, 
it is worth noticing that the BOA correctly describes the internal degrees of freedom in 
terms of translationally-invariant coordinates: q 2 and 23 — {k\ 3 x\ + /c23^2)/(/ci3 + ^23)- 

Finally, we show that we can obtain the BO eigenvalues from the exact solution (1241) . 
Because the BOA is based on the fact that the nuclear masses are much greater than the 
electronic ones, we substitute mi = Ui/X and m 2 = u 2 /A into the roots w% and iy 2 of the 
characteristic polynomial ( |23i) and rearrange the results in order to remove the poles (for 
example, multiplying numerator and denominator by A 2 ). If the undetermined quantities u\ 
and U2 are of the same order of magnitude as m 3 then A < 1. Then we expand both roots 



S 



in Taylor series about A = 0: Wi = w^o + u^iA + . . .. After tedious algebraic manipulation 
of the equations (greatly facilitated by any available computer algebra software) we obtain 

ki3 + k 23 . k\ 3 u 2 + k\ z u x 

W\ — 1 — -A + . . . 

m 3 Wi« 2 (&i3 + k 23 ) 
fcig + k 2 3 kf 3 m 2 + fc^mi 
m 3 mim 2 (fci 3 + k 23 ) 

(Ul + U 2 ) [fci 2 (fcl3 + ^23) + ki 3 k 23 ] . 
W 2 — — r A + . . . 

uiu 2 {k 13 + k 23 ) 

(mi + m 2 ) (k 12 k 13 + k 12 k 23 + k 13 k 23 ) 

+ . . . (35) 



mim 2 (A;i3 + k 23 ) 

The first contribution to w\ gives the large frequency due to the fast motion of the electron 
(second term in the right-hand side of equation (|33|) ). and the dominant contribution to w 2 
provides the frequency for the slow motion of the nuclei (third term in the right-hand side 
of equation (|33|) ). If, in addition, we take into account that rs and q% ~ R when 
we neglect the electron mass, we clearly appreciate that the BOA yields the leading terms 
of the expansion of the exact result in negative powers of the nuclear masses (or the ratio 
of the electron to nuclear mass), even for the contribution of the motion of the center of 
mass. We have decided to express that ratio conveniently as A = Ui/rrii but other forms are 
possible. 

An alternative mathematical strategy for obtaining the expansions of the roots of the 
characteristic polynomial is to substitute w = + w^'X + . . . (and, of course, m, = itj/A) 
into the characteristic polynomial ff23l) and expand the resulting expression in powers of A. 
At order zero we obtain two roots u>j 7^ and = 0, and the coefficients of greater 
powers of A yield further corrections. The result should also be equation ( 1351) . 

V. IDENTICAL NUCLEI 

The particular case of identical nuclei leads to simpler expressions because 

mi = m 2 , k 13 = k 23 (36) 



We easily obtain 



and 



2 _ k 13 (2m 1 + m 3 ) 2 _ 2k 12 + k l3 
m\m 3 mi 

^2m 1 m 3 (2q 3 - q 2 ) y /2m 1 

2/2 = — /0 — , 2/3 = -^—Q2 38 

2^/2m 1 + m 3 2 
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We appreciate that P\ 2 q\ = <?i, P12Q2 = —Q2 and = <?3 — <?2) where P\ 2 is the 



permutation operator that satisfies Pi 2 f(xi,x 2 ) = f(x 2 ,xi). Therefore, Pi 2 y 2 = V2 and 



In this way we can build symmetric and antisymmetric states (including the spin) for boson 
and fermions, respectively. 

At first sight it may seem that when choosing the particle 1 as the coordinate origin we 
are violating the quantum-mechanical principle that states that identical particles are indis- 
tinguishable. However, when we express the resulting wavefunctions in terms of the original 
variables we realize that we can take into account the correct permutational symmetry ex- 
plicitly, and, therefore, there is no violation of that principle. That the identical particles are 
treated exactly in the same way is more clearly seen in the form of the variables that appear 
in the exact square-integrable eigenf unctions: q 2 = x 2 — X\ and 2q 3 — q 2 = 2x 3 — (xi + x 2 )/2. 
We also realize that the Born-Oppenheimer states exhibit exactly the same symmetry as 
follows from the fact that the nuclear factor is a function of q 2 and $ e (^3; x±, x 2 ) actually 
depends on 2x3 — (xi + x 2 )/2. 

The eigenfunctions are somewhat complicated to write them down explicitly here. How- 
ever, the correlation functions for the ground state are not so cumbersome. For example, 



gives us the probability of finding one nucleus at x 2 if the other one is at x\. The fact that it 
exhibits a maximum at x\ = X2 is a consequence of the unrealistic attractive nucleus-nucleus 
interaction of our model. This equation clearly shows that we are correctly treating both 
nuclei as indistinguishable particles. 

If, on the other hand, we integrate over the coordinate of one of the nuclei we have 



Pi2U3 = — 1/3) so that we conclude that 




(39) 




(40) 





\/2mim3UJiU} 2 2m i m 3"i"2 

'' 0<ry-t . . irt -L- <ry-> {/.>-, _L 



(x 1 -x :i ) 2 



(41) 



yV [2m 1 cu 2 + m 3 (cui + uj 2 )\ 
that shows the coupling between the electronic and nuclear motions. 



The BOA yields remarkably similar expressions: p BO (xi,x 2 ) = p(x\,x 2 ) and 



bo/ \ \J2mxm z uj x u) 2 _2_ 

>* u ( Xl ,x 3 )= e 2 " 





(42) 
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Notice that the only difference between the exact (HT1) and BO (1421) correlation functions is 
the neglect of the small frequency with respect to the large one: u\ +u>2 wwi. We appreciate 
that for most purposes the BOA gives a correct description of the system behavior, at least 
according to the simple harmonic model discussed here. 

VI. FURTHER COMMENTS AND CONCLUSIONS 

In this paper we have applied the BOA to an exactly solvable model for a diatomic 
molecule with one electron. We could thus show that the clamped-nuclei approach provides 
the leading terms of the expansion of the eigenfunctions and eigenvalues in powers of the 
ratio of the electron to the nuclear mass. In addition to the internal degrees of freedom we 
also considered the motion of the center of mass that is commonly omitted in most textbooks 
on quantum mechanics and quantum chemistry^. Although molecular physicists are more 
interested in the internal (spectroscopic) degrees of freedom we explicitly considered the 
factor that corresponds to the motion of the center of mass for pedagogical purposes. As 
discussed above for the virial theorem, only the internal degrees of freedom should be taken 
into account for the calculation of expectation values and transition probabilities. The reader 
may find a rigorous discussion of the treatment of the motion of the center of mass for actual 
molecular systems elsewhere^. 

The harmonic potential-energy function of the present model is unrealistic, but such 
a choice allows us to solve the Schrodinger equation analytically and thus compare the 
exact solution with the BO one. Unlike earlier oversimplified pedagogical models used to 
discuss the BOA-^, our harmonic molecule contains the minimum number of particles to be 
a molecule and is therefore more realistic from this point of view. For example, we could 
address the interesting case of identical nuclei. 

The harmonic model is not suitable for the discussion of the Franck-Condon principle 7 
because all the electronic curves exhibit exactly the same nuclear equilibrium distance (q2 = 
0) and frequency as shown by equation (j28|) . Therefore, the value of an overlap integral is 
either one or zero because the harmonic-oscillator eigenfunctions are orthogonal. 

The three-particle harmonic model is also useful for the study of the mass polarization in 
atoms^. In fact , if we choose the particles 2 and 3 to be electrons we have a one-dimensional 
harmonic version of the Helium atom. We can thus, for example, estimate the effect of 
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neglecting the motion of the nucleus, or discuss the isotope effects on the atomic properties. 

Our discussion of the harmonic diatomic molecule may also serve as an introduction to 
the non Born-Oppenheimer calculation of molecular properties^ that requires an adequate 
separation of the motion of the center of mass in the way shown in Sec. HT1 
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